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Abstract 



Electrodynamics of charged scalar bosons and spin 1/2 fermions is studied at non-zero 
temperature, chemical potentials, and possible Bose condensate of the charged scalars. De- 
, bye screening length, plasma frequency, and the photon dispersion relation are calculated. It 

is found that in presence of the condensate the time-time component of the photon polar- 
ization operator in the first order in electric charge squared acquires infrared singular parts 
proportional to inverse powers of the spatial photon momentum k. 

r-| ; 1 Introduction 



Gauge theories in medium with non-zero temperature T and chemical potential /i are of interest in 
^> ' plasma physics, in particular for quark-gluon plasma [HE], astrophysics PI'S] and cosmology [11 [5]. 

. While in heavy ion collisions, where quark-gluon plasma might presumably form, and in stellar 

\ environment chemical potentials are typically large, /U > T, in cosmology it is usually assumed 

' that fj, <^ T. This assertion is based on the tiny value of the observed baryon asymmetry, 

P = {Nb — Nb)/N^ « 6 • 10"^'^ and the bounds on the cosmological lepton asymmetry [6]. 
However, it is not excluded that in the early universe some charge asymmetries, including the 
lepton one, might be large [7]. 

Bearing in mind possible applications to cosmology and, maybe, to other mentioned above 
fields, we consider in this paper a simple case of an Abelian gauge theory, namely electrodynamics 
^ ' and focus on the photon propagation in hot degenerate plasma consisting of charged scalar bosons 

■ with mass ms and fermions with mass mp. This problem is interesting by itself from theoretical 

point of view and may be useful for understanding phenomena in heavy nuclei or neutron and 
quark stars as well as in the early universe. A lot of work has been done in gauge theories at 
non-zero T, see e.g. books [8] and references therein. The usual spinor electrodynamics (SpQED) 
at T 7^ was studied ages ago in classical plasma physics, see e.g. [9l, and later many works were 
dedicated to the investigation of SpQED in the frameworks of relativistic quantum field theory. 
Much less attention was devoted in earlier years to the scalar quantum electrodynamics (SQED) 
in medium by an evident reason that no fundamental charged particles were known. However 
subsequent theoretical discoveries of supersymmetry (SUSY) and spontaneous symmetry breaking 
(SSB), and possible existence of fundamental scalars in some extension of the Minimal Standard 
Model revived interest to SQED at non-zero temperatures, for a review see e.g. [10]. 

Usually SQED was studied at non-zero T but with vanishing chemical potential of bosons, 
because the medium was assumed to be electrically neutral. There are some works that invoke a 
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chemical potential for bosons and also deal with Bose condensation. They are typically interested 
in the relationship between spontaneous symmetry breaking and Bose-Einstein condensation, as 
in im [T2] and do not consider aspects regarding photon propagation in the medium. On another 
hand, the propagation of a heavy quark through a quark-gluon plasma in presence of a quark 
chemical potential is considered in [13j . 

In our work we do not confine ourselves to the restriction /i = and moreover we consider, 
in particular, the maximum allowed value of the chemical potential of scalar bosons which is 
equal to their mass, /ig = ^B, and at which they may Bose condense. The SQED with Bose 
condensate of charged scalar bosons was not considered in the previous literature, except for the 
two recent papers [MKTS], where somewhat different issues were studied. In the above mentioned 
papers and in this one the medium was assumed to be electrically neutral because of the assumed 
compensation of electric charge density of bosons by the opposite charge density of fermions. 
In ref. [14J a conductive spherical object of a finite size, e.g. a star, was considered and the 
neutralization of the interior was achieved by forcing an excessive charge to the surface. In 
ref. [15] as well as here the charge neutrality of the medium was assumed ab ovo i.e. by the initial 
condition that the total (conserved) electric charge densities of bosons and fermions are equal by 
magnitude and opposite by sign. The restriction of the overall charge neutrality is not necessary 
and cosmology with non- vanishing electric charge density might exist (see below). 

A similar situation to that discussed in the present work was studied in ref. [16] where a 
medium with charged scalars and fermions, each with non-zero charge density, was considered. 
The paper was dedicated to different problems but not to the calculation of the polarization 
tensor of photon in such a medium, which is a primary aim of the present work. The charged 
boson state in ref. [16] was described in a different from ours gauge which is not convenient for 
the perturbation theory which is used in this paper. Since the expression for the polarization 
operator is gauge invariant in QED, it should be the same in any gauge. Nevertheless we have 
explicitly verified that our result for the photon polarization operator remains the same in the 
gauge of ref. [16j as well. 

There are two standard approaches to field theory at nonzero temperature and chemical 
potentials: either imaginary time method (Matsubara formalism) or real time method (Schwinger- 
Keldysh formalism). The former is applicable only to the case of systems in thermal equilibrium, 
while the latter is valid for any state of the medium. For a review of these methods, see e.g. 
ref. [17] or books [8]. Below we calculate the photon Green's function for arbitrary medium in a 
physically transparent and simple way taking the expectation values of relevant operators not only 
over vacuum but over any state of the system, either it is a collection of particles with arbitrary 
occupation numbers or a coherent field state (the latter is possible for bosons only). We calculate 
the photon polarization operator in the medium and derive from it the electrostatic (Debye) 
screening length, the plasma frequency and more generally the photon dispersion relation. In all 
previously studied cases our results coincide with the known ones. 

When setting the boson chemical potential equal to its limiting value, ^ib = ^b, we found a 
surprising equation of motion for the photon: the correction to the time-time component of the 
photon polarization operator at zero frequency contains in addition to the free term Gqq = k"^ 
and the well known momentum independent term, which is inversely proportional to the Debye 
screening length squared if < m^, two additional terms, one of which appears when the 
bosonic chemical potential reaches its maximum value, ^ = niB, and is inversely proportional 
to the photon momentum, i.e. 1/k, and another term which appears when bosons condense and 
which is inversely proportional to the photon momentum squared, i.e. 1/fe^. In other words the 
Maxwell equation for the time component of the potential at zero frequency (stationary case) 
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and in the low momentum limit is modified as: 

- [k^ + e^{ml + ml/k + ml/k'^)]AQ = Jq, (1) 

where rrij are some constant parameters expressed through the characteristics of the medium 
and Jo is the electric charge density of the medium, see below eq. (f2T]) . As is mentioned above, 
we assumed that Jq = to avoid "infrared" problems, though this restriction is not necessary, 
especially in cosmology. The 1 //c^ term was also found by G. Gabadadze (private communication, 
work in progress). We thank him for sending the draft of his work prior to publication. In the 
paper by Gabadadze a zero temperature case is considered and because of that 1/k term does 
not appear. As shown in the present work, the 1/k term leads to a power law decrease of the 
screened potential, eq. ([59|) . 

To study the screening phenomenon we take, as usually, Jq to be a point-like test charge, Jq ~ 
g(5(3)(k). Evidently the Debye screening mass in this case is not equal, as usually, to erriQ but to 
more complicated expression involving mi and m2 and possessing a non-zero imaginary part which 
leads to an oscillating behavior of screening superimposed on the exponential decay, see eq. (I58|) 
. Similar oscillating behavior of the screening was also observed in ref. [15]. The oscillations of 
the screened electrostatic potential was found long ago in non-relativistic fermionic plasma, see 
e.g. book [18]. They are known as Priedel oscillations. More recently a similar phenomenon was 
investigated in ref. [19j for relativistic QED and QCD0. The oscillating behavior due to fermions 
appears only in higher order in the coupling constant, a = e^/An, while oscillations due to Bose 
condensate appear in the lowest order in the coupling. 

The paper is organized as follows: in Sec. [2] we introduce the photon polarization tensor 11'^'^ 
and the techniques we used to work it out. In Sec [3] we use the appropriate components of 11'^'^ 
to calculate the Debye screening mass rriD and the plasma frequency uip. Some limiting cases, 
which do not include condensate, are analyzed to check compatibility with previous results found 
in literature. In Sec. [Hwe calculate the Debye mass in an explicitely fixed gauge. Finally in Sec. 
[5] our conclusions are presented. 

2 Perturbative calculations of the photon polarization operator 

The Lagrangian of interacting electromagnetic field and charged scalar and fermion fields with 
masses and mp respectively and with opposite electric charges ite has the form: 

£ = - m||,/.|2 + 1(9^ + i eA^)(t>\^ + - e4 - mp)^. (2) 

The Lagrangian is symmetric under the gauge transformations: 

(j){x) ^ eyiY>[iea{x)](j){x), iIj{x) ^ e^Y>[iel3{x)]il){x), ^ - 9^(a -h /?), (3) 

This implies an existence of 2 conserved currents and charges, which we can choose as the scalar 
and the fermion number. The Lagrangian ([2]) leads to the following equations of motion for the 
involved fields: 

— m)ip{x) = e4'(p{x) (4) 

{d^df' + m^)<P{x) = J^x) (5) 
*We thank G. Gabadadze for indicating these references. 
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(6) 



where the currents J are defined as 



J4,(x) 



-I e 



{x)^^A^'{x)A^(x)^(x) 



(7) 



J^ix) 



-I e 



{<\>\x)d^4>{x)) - (d^<\>\x))4>(x) 



+ 2e^A^{x)\4)(x)\^ - e^(x)-i^^{x). (8) 



Here Fp,y = d^A^ — dyAp, and Q is the total electromagnetic current of bosons and fermions 
in the coordinate space. 

The key quantity which determines the photon propagation in plasma is the photon polariza- 
tion tensor which we will calculate perturbatively. When doing this kind of calculations for 
massless fields, infrared singularities may arise and to regularize them one should use the resum- 
mation techniques - see e.g. [8] and references therein. Nevertheless it is safe to use the standard 
perturbative solution when the scalar and fermion field masses are not negligible as in the case 
we are considering. Moreover, the infrared singularities in Abelian theories are much milder than 
those in non- Abelian ones where the correspondence between the order of the perturbative series 
and the power of the coupling constant e is lost, see e.g. discussion in ref. [H |20]. Here we 
consider only Abelian QED and since it is not infrared dangerous, we neglect the resummation. 

We formally solve operator equations (jl]) and ([5]) as: 



V'(x) = iJo{x) + J d'^yGpix - y)e4{y)ij{y) 
where the zeroth order fields satisfies the free equations of motion: 

{^^,^^' + m%)Mx) = 0, (i^ - mF)Mx) = 

and are quantized in the usual way 



(9) 



(10) 



il^oix) 



d^q Imp 



y(2^V E 



a(q) exp '^^ +6'!" (qj exp 

Cr(q)ur(q) exp"*^"^' +4(q)ur(q) exp*"?^ 



fill 



In eq. (jlip a^'^\ b^'^\ c^^^ and d^^^ are the annihilation (creation) operators for scalar and spinor 
particles and antiparticles. The Green functions in eqs. ([9]) are the usual Feynman Green 
functions having the integral representation: 



where 



fc2 



+ ie 



, and Gpik) = 



fc2 



mir + le 



(12) 



(13) 
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Now we can substitute eqs. ([9]) into eq. ([6]) with the currents given by eqs. ([7]) and ([8|). For the 
calculations up to the second order in the couphng constant e, i.e. up to e^, it is sufficient to 
include into Jtf, in eq. ([7|) only terms of the ffist order in e, that is 



As a result we obtain 
a^F^^(x) = -ie 



+i e 



^^A^'ix) + 2A^ix)^^' 



(14) 



d^yGsix - y)J^o{y) 

d^yGsix - y)J<j,„{y) 
- eipo{x)-f^ J d'^y Gf{x - y)e4{y)7p{y) 
+ 2e'A^^{x)\Mx)f- 



t 



e'ipo{x)j^''ifjo{x) 

d^yGsix - y)J4>oiy) 
d'^yGsix - y)J4>o{y) 

d^yMy)4iy)GUx-y) 



t 

Mx) 
i^i^oix) 

(15) 



To derive the Maxweh equations with the account of the impact of medium on the photon 
propagator we have to average operators (j) and V' over the medium. The ffist term in eq. (jlSp . 
Unear in e, is non-zero if the medium is either electrically charged or possesses an electric current. 
The products of creation-annihilation operators averaged over the medium have the standard 
form: 



(at(q)a(qO) 
(ct(q)c(qO) 



/B(i?,)<5(3)(q_q'), 
/^(Sp)5(3)(q 



(a(q)at(q')) = [l + /B(i?p)]<5(3)(q_q'), 
(c(q)ct(q')) = [l-/^(Sp)]5(3)(q-q'), (16) 



where fF,B{Eq) is the energy dependent fermion/boson distribution function, which may be 
arbitrary since we assumed only that the medium is homogeneous and isotropic. We also assumed, 
as it is usually done, that non-diagonal matrix elements of creation-annihilation operators vanish 
on the average due to decoherence. For the vacuum case fF,B{E) = and we obtain the usual 
vacuum average values of aa) and a^a, which from now on will be neglected because we are 
interested only in the matter effects. As a result we obtain linear but non-local equation for 
electromagnetic field A^{x). It is convenient to perform the Fourier transform 



d^x 
J2nf 



■ exp 



-ikx 



Ai'ix). 



Finally we find that the field A^(k) satisfies the equation 

[kf'kpg'"' - k^'k'' + W'^ik)] A,{k) = J^'{k), 



(17) 



(18) 



which is equivalent to the photon equation of motion ([6]) but in momentum space. 

Thus the photon polarization tensor which contains contributions from charged bosons and 
fermions, Il^,^{k) = Ilj^^{k) + I[j^j^{k), and the electromagnetic current J'^ involved in eq. ([TH]) 
are explicitly found in the lowest order: 



d\ 



[fB{E)+fB{E)] 



I {2q-k)^{2q-k)^ , I {2q + k) ^{2q + k) 



(q - kf 



+ 



m 



B 



{q + kf 



m 



B 



,(19) 
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Qvik + q)fj. - qPkpQ^u + q^ik + q)i 



{k + qf 



+ 



Qiy{q - k)p + q^kpQ^y + Qfiiq - k)y 



ml 



{k-qf 



ml 



, (20) 



Ju 



(2^ 



■ exp 



d^q 



9m 



(2^)3 E 



fB{E) - fB{E) - 2[fF{E) - fF{E) 



(21) 



In eqs. ()19p - (j2ip k^^ = (a;,k) and q^ = (i^', q) are respectively the photon and the scalar /spinor 

four momenta, fq and fq are the particle (antiparticle) distribution functions and g^u = (H )• 

We assume the following charge convention: the bosons have electric charge +e, while fermions 
have electric charge — e. Of course the charge of antiparticles has the opposite sign. It is worth 
to stress that the total 11'^'' as well as its bosonic and fermionic components separately satisfy the 
transversality condition k^Yl^^y = 0. The last term in eq. ()19p describes the contribution from 
the tadpole diagram and coincides with that found in ref. [21J. Evidently without this term the 
tranversality condition would be violated. 

In what follows we use 11'^'^ to derive plasma frequency, Debye screening mass, and dispersion 
relation for the photon propagation in plasma. 



3 Photon propagation in plasma 

An electric field of a test particle is known to be screened by the plasma polarization. In the 
static case i.e. for lij = it follows from eq. (jlSp : 

(|k|2-noo)Ao = -g, (22) 

where g is a (small) charge of the test particle. When IIoo is independent on the photon momentum 
k, eq. (j22p becomes the usual equation for a scalar field with mass m = -y/— IIoo- In this case 
which is usually realized at least for small |k| the Debye mass mz) is equal to a/— IIoo and the 
electrostatic potential turns from the Coulomb to the Yukawa one: 

q exp(-mz)r) 

Hence the Debye mass coincides with the inverse of the screening length. More generally the 
Debye mass is defined as the position of the poles of the inverse operator (|kp — IIoo)"^ and in 
presence of the Bose condensate the simple relation mo = \/—^oo{^ = 0) becomes invalid (see 
below) . 

There may be some ambiguities in definition of the screening potential at higher orders in 
electric charge, and higher, and in non-Abelian gauge theories where the result may violate 
gauge invariance, as discussed in ref. |22j . However, we consider here only Abelian U{1) theory 
and in the lowest, order. Moreover, as one can see below, we have calculated the asymptotics of 
the screening potential at large distance in a self-consistent way by the location of the singularities 
in the complex A: -plane, which are situated at non-zero k. 

Another quantity of interest is the plasma frequency iVp, which enters the dispersion relation 
of electromagnetic waves propagating in plasma. It is defined as the limit of k — > of certain 
space-space components Hy. Using the transversality condition, k^Ilf^u = 0, we can decompose 



6 



the photon polarization operator in a medium in terms of two scalar functions, a{uj, |k|) and 
6(a;,|k|): 

% = a(5.,-Ml)+,^, = ^5, noo = ^6. (24) 

We have assumed here that the medium is isotropic and so a and b depend only on the absolute 
value of the photon momentum. 

One can check from the presented above explicit expressions (|19p and (|2U|) that a = b m. the 
limit k = 0. In other words, in this limit liij ~ 8ij. The plasma frequency is determined by the 
equation: 

= 6(cj, |k = 0|) = a(cj, |k = 0|). (25) 

Since the distribution functions depend only upon the energy, the integral over angles in eqs. 
(jl9p and (j20p can be taken. In particular for time-time component of the polarization tensor in 
the limit of u; = we find: 



e 



2 roo A^^'i 




Here and in what follows k and q are respectively the absolute values of the spatial component 
of the photon and charged particle momenta. The argument of the logarithm comes in absolute 
value because the Green's functions in the perturbative expansion appear in the combinations 
G{q + k) + G*{q — k). So the imaginary part of IIoo in the considered limit vanishes. 

It is even simpler to calculate the space-space components, Iljj in the limit of zero photon 
momentum. A: = 0: 

= .J ,., ]{fB + fB), (28) 




nfj = / \'; {fF+fF). (29) 

It is clear from these expressions that in the limit A; = functions a and b (|24p are equal. 

At u; = 2m the polarization operator acquires a non-zero imaginary part which corresponds to 
the threshold of two charged particles production by the photon. For massless charged particles 
the threshold is at w = and one has to take into account that the effective photon mass in 
plasma is non-zero (it is essentially the plasma frequency). This can be done using resummation 
technique. Accordingly the position of the threshold moves a little, as ~ eT. This is not essential 
for our consideration. Moreover electrodynamics with massless charged particles has serious 
infrared problems. 

Let us now apply these results for the calculation of the plasma frequency and the Debye 
mass in some special cases which have been considered in the literature. 

Massless SpQED and SQED with vanishing chemical potentials and high temperature, T ^ 
mB,F- In thermal equilibrium the distributions of bosons and fermions and their antiparticles 
with zero chemical potentials have the usual Bose-Einstein or Fermi-Dirac form: 

fB{E,T) = fB{E,T) ^ 



exp(^/r) - 1' 
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ME,T) = ME.T) = ^^^^^^^^^ , (30) 

where in high temperature hmit we can neglect the particle mass i.e. we can assume -E = |q|. In 
this special case IIoo does not depend on the photon momentum k and so = \/—IIoq. The 
integrals in eqs. (jl9p and (j20p can be easily taken and we find for the contributions from bosons 
and fermions respectively: 

"T'Db("^b = 0) = m^pimp = 0) = ^e'^T'^, 

w|,b("ib = 0) = ujlp{mF = 0) = le^T^. (31) 

y 

Hence the total Debye mass and the plasma frequency for the Lagrangian ^ are 

mlimB = = 0) = ^ e^T^ (32) 

LOp{mB = mp = 0) = n (33) 
9 

These expressions coincide with the published results in the lowest order in the electromagnetic 
coupling, e^, see e.g. [TOl [8]. 

In the case of relativistic fermions with non-zero chemical potential fi we obtain: 



mUmp = 0) = — + ^ . (34) 



r2 

h — 

This is the result found in ref. [23]. 

Massive SpQED with non-zero chemical potential and low temperature, T ^ 0. In thermal 
equilibrium the distributions of fermions and their antiparticles are in this case the Boltzmann 
distributions: 

fp{E, /X, T) = e(^-^y^, fp{E, T) = e-(^+^)/^ (35) 

and once again IIoo does not depend on k, so it coincides with the Debye mass squared. Hence 
in the case of relatively small chemical potentials, fi < m, we obtain 

noo = ^ (^) e--/^ cosh(^/r). (36) 

In the case of strongly degenerate, /i > m, nonrelativistic fermionic plasma, the contribution of 
anti-fermions may be neglected, while the fermion distribution function has the form: 

The chemical potential can be expressed through the number density of the fermions: 



exp(/i/r) 



TT 



2 -q^/2mFT 



up = ' I dqq'e-1 r^""^' . (38) 



There is a factor 2 in the above expression which counts two spin degrees of freedom. 
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Correspondingly the Debye screening mass for nonrelativistic fermions is 

mo = -j^, (39) 

which coincides with the classical result, see e.g. book [9]. 

Analogously we find from eq. ()29p the plasma frequency for nonrelativistic fermions: 

ujI = , 40 

which is also the classical result. 

We have done these simple exercises to check the validity of our results for li^u comparing it 
to the known cases. Now we will turn to the calculations of the photon polarization tensor in the 
medium with Bose condensate of charged scalars, which to our knowledge was never done before. 

3.1 Bose-Einstein condensate 

If the density of particles is not equal to the density of antiparticles, the equilibrium distributions 
are modified by introduction of chemical potentials. The distribution functions take the form: 

feAE. T) = (41) 

where the plus sign stands for fermions and minus is for bosons. In equilibrium the chemical 
potentials of particles and antiparticles have the opposite signs, p. = — /i. Evidently chemical 
potential of bosons is bounded from above, < f^B- 

If the asymmetry between bosons and anti-bosons reaches so large value that jiB = 
is not sufficient to ensure such a large asymmetry, Bose condensate would be formed and the 
distribution functions take the form: 

f^S\E,C,T) = ^ +C6{cD^fB{E,mB,T) + C6{cD, (42) 

exp[(^ - mB)/T\ - 1 

fB{E,-mB,T) = \ (43) 

exp[(£; + mB)/T\ - 1 

where a new constant parameter C is the amplitude of the condensate. 

The distributions in eq. (j42p are the stationary solutions of the kinetic equation if and only if 
fJ,B = iTT-B (or fj,B = itt-b)- Distribution (j42p with C ^ describes non- vanishing number density 
of bosons in a vanishingly small momentum interval near q = 0. In this sense it is a classical 
field configuration (j){t) = (/)o exp(imBt), where the field amplitude is related to the condensate 
amplitude as 2(/)gm = C/{2tt)^. On the other hand, it evidently can be interpreted as a collection 
of particles at rest. 

Once the boson charge density J^i^ is fixed, it is possible to calculate the critical temperature 
Tq when the phase transition takes place and the Bose condensate is formed. This temperature 
can be calculated from eq. (f2T]) imposing 

Jcf {^,B = mB,C = 0) = (Tc), (44) 

from which follows the equation 

Tc = (45) 
V eniB 
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which coincides with what we have found in the hterature - see e.g. ref. |TT] . 

We can express the amphtude of the condensate through the electric charge density of the 
plasma. We impose the condition that the total electric charge density of fermions and bosons is 
zero. The condition of vanishing electric charge density may be not fulfilled e.g. in cosmology [24j, 
but we postpone this case for future consideration. So we assume that = in eq. (pT]) . 
Expressed through the particle occupation numbers /s,_f this condition reads: 



d^q [fB -Ib- 2{fF - 2fF)] 



0. 



(46) 



Here we skipped the arguments of fF,B, which are explicitly given in eqs. (|4m43p . The factor 2 
in front of the fermionic distribution functions is related to two spin degrees of freedom. 
Hence the amplitude of the Bose condensate C for globally neutral plasma is equal to: 

C = -iTT f dqq^ [fBiEB,mB,T)-fBiEB,-mB,T)-2fF{EF,tiF,T) + 2fF{EF,tiF,T)]. (47) 



Here Eb,f = y F + 1 • "^^^ magnitude of the fermionic chemical potential, fiF, is determined 
by the charge density of fermions. In the model under consideration it is a free parameter. We 
can imagine another more realistic model, when the bosons have the Yukawa coupling to electrons 
and neutrinos, This coupling leads to the process (/>"'' <-> + e^. If the lepton number 

density is sufficiently large, chemical potential of (j) may reach the limiting value, fj,B = w-s, and 
(/)-bosons would condense. 

Now we can use eqs. p^ . (p6|) and (p7|) to calculate Iloo{uJ = 0, k) which is needed to work 
out the electrostatic potential of a test charge in the plasma with Bose condensate: 



C 



1 + 



4m^ 



fc2 



B 



e 

Hoo (w = 0,k) = ——-r 

{2TTy rriB 

+ ^ ^^^^ [fB{EB,mB,T) + fBiEB,-mB,T) 
+ £^ 1^°^ ^[fF{EF,^^F,T) + fF{EF,-l^F,T)] 



1 + ^ln 

kq 



2q + k 



2 + 



2q-k 

{AEl -k^) ^^ 2q + k 
2kq 2q-k 



.(48) 



It is rather surprising that the first term which appears only in the presence of the Bose 
condensate has a singularity, l/fc^. However this is not the only singular in k term. There is also 
a singularity 1/k arising from the bosonic (but not anti-bosonic) contribution to Hoo if = 
even in the absence of condensate, i.e. for C = 0, which appears from the integration region 
where q is smaller or comparable to k. For small q the distribution function is infrared singular. 



fB{EB,mB,T) w 2mBT/q^ 



(49) 



Usually this singularity is not dangerous because it is canceled by the integration measure, ~ q^. 
However, there is another factor to worry about. The logarithmic term behaves as k/q for q > k 
and as q/k for q < k. So the integral is finite but has 1/k singularity. Indeed we can separate 
the integral into two parts < q < k/2 and k/2 < q < oo. It is convenient to introduce for 
the first part the new integration variable x = 2q/k, so < x < 1. For the second part we 
introduce y = k/2q, so y runs in the same limits, < y < 1. In the limit of small k we can 
expand Eb ^ tub + /c^x^/Sm^. Correspondingly: 



exp 



Eb - ruB 



1 « k'^x^/SruBT. 



(50) 
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So the integral is reduced to 



^ dx ^ f 1 + X 

— In 

X \1 — X 



vr 

T 



(51) 



The same contribution comes from the second part of the integral. So finally we obtain for the 
singular in k part of the photon polarization tensor in the case of /u = and C = 0: 



n 



00 



e^m%T 
2k 



(52) 



Numerical calculations without expansion of the energy and the exponent gives a very close result. 

Thus for zero frequency and in the limit of small k the time-time component of the polarization 
tensor has the form: 



n, 



00 



2T^ 



m'lT 
+ -5- + 



1 



C 



2k (27r)3 ruB 



1 + 



Am 



B 



k^ 



(53) 



This form of IIoo determines the position of the singularity points in the complex k-plane at low 
k, which in turn governs the asymptotics of the screening potential at large r. For simplicity 
we took here the high temperature case, T ^ mp^B when the k- independent terms have a very 
simple form but the analysis can be done for any T. The general form of IIoo is presented in eq. 
(dl). According to the result above: 



mn 



mo 



2rV3 + C/[{2TTfmB] 

m|r/2 

4CmB/(27r)3. 



(54) 



To determine the Debye screened potential we have to make the Fourier transformation of 
the inverse of the photon equation of motion k"^ — Uqq = 0: 



Uir) 



dkk'^ 



sin kr 



d^k exp(ikr) q 
(27r)3 fc2 _noo(A;) " 2^2 j,2 _ n^^^j.^ 



(55) 



If IIoo is an even function of k, as it is usually the case, the integration over k can be extended 
to the interval from — oo to +oo and the integral can be taken as a sum over residues of the 
integrand. For example, if the term proportional to 1/A; can be neglected (low temperature case), 
IIoo is evidently even and its poles can be easily found: 



kj 



e^ml 



\ 



4^4 



e^m, 







9 4 



±iem2- 



(56) 



The last approximate equality is formally true in the limit of small e. 

Thus in presence of Bose condensate the "Debye" poles have non-zero real part: 

kj = ±\/em2 exp(ibi7r/4) = kj + ik" , (57) 

while A:' = in the usual case. Non-zero k' leads to the oscillating behaviour of the potential 

exp(— ^e/2m2r) cos{y^e/2m2r) 



U{r), 



(58) 
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If the term proportional to \/k is present in IIoo the calculations of the potential are slightly 
more complicated. Now the integration path in the complex k plane cannot be extended to — oo 
but the integration should be done along the real axis from to oo, then along infinitely large 
quarter-circle, and along the imaginary axis from — oo to 0. The result would contain the usual 
contribution from the poles in complex /c-plane and the integral over imaginary axis. The former 
gives the usual exponentially decreasing potential, while the latter gives a power law decrease: 

U qm\/{e^mlr^). (59) 

Notice that the potential is inversely proportional to the electric charge squared. 

These strange results possibly indicate that the approximation of the uniform density plasma 
is not accurate near the test charge in the presence of Bose condensate and the polarization of 
the medium is very strong, which is rather natural for particles at rest. 

It is worth to notice that Debye screening is non-analytic in e, i.e. ~ e^/^. 

On the other hand, there is no anomalous behavior for the propagating modes, i.e. for Iljj. 
In particular the plasma frequency in the presence of the condensate acquires simply a constant 
contribution 5uj^ = e^C j {2'Kf'mB- 



4 Calculations in another gauge 

Though our result is gauge invariant and for the calculation of the Debye screening length we do 
not need to fix the gauge, still it may be instructive to make the calculations in the gauge used in 
ref . [16] . The homogeneous state of charged scalars with non-zero charge density was described 
in this paper as 

(j) = (pQ = const, (60) 
<) = (mB/e)<5°. (61) 

It can be easily seen that this solution is a gauge transform of the zeroth order state (in coupling 
e) used in the present paper: 

(j) = (j)oexp{imBt), (62) 
Af, = 0. (63) 

Evidently such a state of (p describes a collection of bosons at rest i.e. of a Bose condensate. The 
electric charge density of the condensate can be read off equation ([8]) and is equal to: 

jP = 2emB\M^. (64) 

Expressions (|60|6ip and (|62|63p lead to the same result for the electric charge density. Comparing 
it with the charge density described by the equilibrium distribution (fl2]) we find that we have to 
identify: 

C7/(27r)3 = 2mB\M^. (65) 

Perturbation theory is less convenient in gauge (j60l61|) because of large value of the back- 
ground potential Aq'^ ~ 1/e. We need to make the expansion: 

^M-^-^M + ^M' (66) 

= (/>0 + (pq, (67) 
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where is the potential of the physical electromagnetic field in the plasma. It may describe e.g. 
Coulomb-like field around test charge when we discuss Debye screening or propagating waves in 
plasma when we talk about plasma frequency. This is equal to considered above when we 
worked in the gauge defined by eqs. (I62|63p . The quantum deviation from the condensed state 
of the scalar field, (j)q, is supposed to be zero on the average. Moreover, we assume for simplicity 
that the plasma temperature is zero and thus ((/>^) vanishes as well (if the vacuum quantum 
fluctuations are subtracted). So (p enters only into description of virtual particles through the 
Green's function. 

The equation of motion of (j)q has the form: 

(d^ - 2imBdo)cj)g = 2emBAocl)o + e^iA^fi^Po + <l)g) + ie [2^^^^ + (d^A'')] (68) 

Only the first term in the r.h.s. of this equation will be essential in what follows. 
It is straightforward to write down the equation for A^: 

d^A^ - [d'A,) = 4^'> - 2emB5l - 2(t>^Re(t>q] + ... . (69) 

Here we retained only terms which are essential for the calculation of the condensate impact on 
the photon propagation in plasma. The missing terms can be easily found from equation (jlSp . 
The first term in eq. (I69|) is the fermionic current. Together with the second term they make the 
total current in the plasma, which we assume as above to be zero, jP - 2emB5l(t)l = 0. The 
remaining term in the r.h.s. can be found by perturbative solution of eq. (I68p : 

(t>q = 2emB J G{x - y) [</.o^o(y) + ••■] • (70) 

As above, we make the Fourier transformation (jl7p and obtain the correction to the time-time 
component of the polarization operator from the field (f)Q: 

OToo = e^ml<l)l/\^. (71) 

Keeping in mind identification (j65l) we find that it is exactly the same result which we have 
found above working in terms of equilibrium distribution with Bose condensate. This finalizes the 
argument that both descriptions (in both gauges) are equivalent and both states (p = (pQ exp(imt) 
and Aq = and (p = (pQ and Aq = m^/e describe the same collection of particles at rest. 



5 Conclusion 

We have calculated the polarization operator of photon in the plasma with charged scalar bosons 
and spin one half fermions. The calculations have been done in the lowest order in electromagnetic 
coupling constant, e. Our technique is slightly different from the standard one, that is from either 
imaginary time method applicable to thermal equilibrium case or real time method which leads 
to the doubling of states. We simply solved perturbatively the operator equations of motion for 
the charged fields The solution in the lowest order typically has the form 

<^i{x) = J dyG{x -y) J i^o), (72) 

where $o,i are the charged fields in zeroth and first approximation in e respectively and is 
the current entering the equation of motion. This solution was substituted into the Maxwell 
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equation for the quantum electromagnetic field and the average of quantum operators ^>o was 
taken over the medium. In this way the effect of medium for arbitrary occupation numbers of the 
charged particles (not necessary equilibrium) can be taken into account. We checked our results 
comparing them with the known cases of Debye screening and the plasma frequency at low and 
high T and nonzero chemical potentials. 

The new part of the work consists of calculation of the polarization operator in the presence 
of Bose condensate of charged scalar particles. Though it is technically simple, there are no such 
calculations in the literature, as far as we know. Physically such a case can be realized if there 
is a significant asymmetry between charged fermions (electrons and positrons). Assuming that 
the global charge of the plasma is zero (though it is not necessarily so) we can find the number 
density of charged scalars in the condensate. To this end the charge asymmetry between fermions 
must be sufficiently large so that the maximum allowed chemical potential of bosons, = is 
not enough to secure the vanishing of total charge. The necessary neutralization can be achieved 
by the charge density of the condensate, aving the distribution fc = C5{q). The corresponding 
correction to the polarization operator has a singularity, i.e. a pole at zero three-momentum of 
photon, 4Ce^m^/(27r)^/c^ + e^m^r/2fc. Probably the origin of this singularity is a large mobility 
of particles with zero momentum under the influence of an external electric field. Energetic 
particles are much less "eager" to screen the test charge. 

Due to these infrared singular terms the Debye screening length becomes parametrically 
shorter, A^) ~ 1/y/e (if C 7^ 0), instead of the standard one Xd ~ 1/e which is true for the 
plasma of charged fermions and/or bosons with < i^T'B- It is noteworthy that the Debye 
screening is now a non- analytic function of the electromagnetic coupling e. The screened electro- 
static potential in the presence of Bose condensate of charged scalars has an oscillating behavior 
superimposed on the exponential decrease, eq. (I58p . The presence of 1/k term which appears at 
non-zero temperature when ^ = niB, leads to weaker, power law, decrease of the screened poten- 
tial, which is inversely proportional to the electromagnetic coupling constant squared, eq. ()59p . 

There are no known charged stable scalar fields, so the considered situation may be unrealistic. 
However, it is not excluded that an extension of the minimal standard model will demand stable 
charged Bose fields. If this is the case, the calculated photon dispersion relation may be of 
interest in cosmology or in dense stellar environment. For example, we may expect condensation 
of di-quarks in quark stars. On the other hand, stability of scalars may be irrelevant because 
even unstable particles may condense in a dynamical equilibrium state. 

Bose condensation of charged scalars can possibly be realized also in the following situation. 
Let there be high temperature plasma of e^, and i> with a considerable excess of e~ over 
and of v over P, i.e. there are large lepton and electric charge asymmetries. We assume that there 
are also charged pions (or some other scalar particles) in the plasma to ensure zero total electric 
charge. The condensed scalars may be stabilized by large chemical potentials of neutrinos, such 
that the decay tt"'" — > e'^v is not allowed because the Fermi states with E = are already 

occupied. Such a state might possibly exist in exotic stars. 
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